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Recent  geomagnetic  measurements  have  shown  that  changes  in 
the  geomagnetic  field  have  deviated  from  th*e  predictions 
of  the  original  International  Geomagnetic  Reference  Field 
(IGRF)  coefficients  and  associated  time  derivatives. 

Accordingly,  a new  set  of  time  derivatives  have  been  de- 
termined with  tlie  suggestion  that  these  be  applied  for 
time  periods  after  1975.  To  ascertain  the  effects  of  the 
changes  of  geomagnetic  field  time  derivatives  on  cosmic 
ray  analyses,  vertical  cutoff  rigidities  were  calculated 
for  selected  locations  using  the  IGRf  1965  field  coeffi- 
cients and  associated  time  derivatives  and  the  IGRF  1975 
field  coefficients  and  the  new  time  derivatives  making 
each  set  of  calculations,  appropriate  for  a 1980  Epoch. 

A comparison  of  these  two  sets  of  calculations  shows  that 
the  uncertainty  in  predicted  magnetic  fields  does  not 
appear  to  be  a serious  limiting  factor  in  the  use  of  cal- 
culated vertical  cutoff  rigidities  for  the  analyses  of 
cosmic  radiation  data. 

1.  INTRODUCTION 

Many  analyses  of  cosmic  radiation  data  are  dependent  upon  the  cutoff  rigid- 
ity of  the  locations  where  data  are  acquired.  For  the  past  fifteen  years 
these  cutoff  rigidities  have  been  determined  by  the  computational  technique  of 
numerically  tracing  the  orbits  of  cosmic  rays  through  a mathematical  model  of  | 
the  magnetic  field.  Although  this  technique  is  considered  as  the  most  reliable 
method  of  determining  cutoff  rigidities,  it  is  also  recognized  that  the  cutoff 
rigidities  determined  in  this  manner  can  only  be  as  accurate  as  the  magnetic 
field  description  utilized  in  the  trajectory-tracing  process.  , 

The  International  Geomagnetic  Reference  Field  (IGRF)  adopted  by  lAGA  in  i 
1969  (lACA  Conmission  Z Working  Group  4,  1969)  is  a composite  of  several  geo-  [ 
magnetic  field  models.  Although  the  model  itself  was  normalized  to  a 1965 
Epoch,  tne  time  derivatives  were  derived  from  actual  field  measurements 
through  1967.  Thus,  the  utilization  of  the  time  derivatives  for  years  prior 
to  l9b5  resulted  in  fairly  representative  models  of  the  actual  field  as  these 
time  derivatives  were  determined  from  magnetic  field  measurements;  however, 
extension  of  the  time  derivatives  into  the  future  (i.e,  beyond  1968)  produced  j 
a "predicted"  field  based  upon  the  assumption  that  the  time  derivatives  them-  i 
selves  were  constant.  I 

Subsequent  measurements  of  the  geomagnetic  field  revealed  changes  that  were 
not  consistent  with  the  predictions  of  the  original  IGRF  time  derivatives.  In 
general,  the  geomagnetic  field  has  been  changing  more  rapidly  than  was  predict* 
ed  in  1968  (Cain,  1975). 

After  consideration  of  the  problems  associated  with  deriving  updated 
models  and  also  considering  the  various  uses  ot  the  models,  the  lAGA  Division 
1 Study  Group  (1976)  recommended  the  following:  that  the  IGRF  1965  raoddl  with 
its  time  derivatives  ^.iiould  be  used  for  the  period  1955-1975,  but  that  a new 
set  oi  time  derivatives  to  give  a reasonable  estimate  of  the  true  rate  of  the 
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I variation  at  Epoch  1975  should  be  utilized  lor  predicting  magnetic 

' jield  coefficients  for  epochs  later  than  1975,  Thus,  a discontinuity  In  the 
j specification  of  the  geomagnetic  field  would  occur  at  1975.  In  this 

• •^»pcr  we  have  Investigated  If  this  discontinuity  In  the  geomagnetic  field 
' ipeelf  feat  Ion , results  in  a similar  discontinuity  in  the  calculation  of 
! yrrclco^  cutoff  rigidities. 


2.  METH(» 

Vertical  cutoff  rigidities  were  calculated  by  the  trajectory-tracing 
technique  using  two  mathematical  representations  of  the  internal  geomagnetic 
((eld  adjusted  to  a 1980  Epoch  as  follows:  (A)  the  IGRF  coefficients  for  Epoch 
1965  (lAGA  Commission  2 Working  Group  4,  1969)  with  their  original  time 
derivatives  applied  for  a span  of  15  years  to  Epoch  1980,  and  (B)  the  IGRF 
coefficients  for  Epoch  1975  (lAGA  Division  1 Study  Group,  1976)  with  the  new 
time  derivatives  applied  for  a span  of  5 years  to  Epoch  1980.  These  two 
representations  will  be  called  the  original  and  adjusted  time  derivatives, 
respectively,  for  the  remainder  of  this  paper,  1 

These  vertical  cutoff  rigidities  were  determined  by  the  method  described 
by  Shea  et  al.  (1965)  using  0,01  GV  Intervals  throughout  the  penurabral  region. 

■These  values  were  obtained  for  the  following  locations;  (A)  cosmic  ray  stations 
In  the  Western  Hemisphere  where  the  vertical  cutoff  rigidity  decreases  with 
time  (Shea,  1971),  and  (B)  selected  locations  on  the  world  grid  of  vertical 
cutoff  rigidities  (Shea  and  Smart,  1975)  where  the  vertical  cutoff  rigidity 
(1)  Increased,  (2)  decreased,  or  (3)  remained  essentially  constant  with 
respect  to  the  1965-1975  interval. 

Once  the  vertical  cutoff  rigidities  were  determined  for  these  locations, 
the  results  were  evaluated  in  the  following  manner: 


TABLE  1 

VERTICAL  CUTOFF  RIGIDITIES  FOR  SELECTED  COSMIC  RAY  STATIONS  IN  LATIN  AMERICA 


Vertical  Cutoff  Rigidity  (GV) 


Station 

Rigidity 

Parameter 

Epoch 

1955 

Epoch 

1965 

Epoch 

1970 

Epoch 

1975 

Epoch 

1980 

Epoch 

1980A 

Mexico  City 

P(m) 

10.06 

9.57 

9.51 

9.44 

9.35 

9 ‘.33 

1 

P(c) 

9.46 

9.12 

8.99 

8.88 

8.70 

8.71 

P(a) 

8.95 

8.77 

8.67 

8.57 

8.47 

8.46 

Buenos  Aires 

P(m) 

10.89 

10.59 

10.44 

10.31 

10.15 

10.13 

P(c) 

10.58 

10.22 

10.03 

9.88 

9.73 

9.69 

P(a) 

10.17 

9.80 

9.62 

9.46 

9.30 

9.28 

Ushuala 

P(m) 

5.80 

5.60 

5.50 

5.42 

5.32 

5.29 

P(c) 

5.67 

5.51 

5.33 

5.29 

5.17 

5.14 

P<8) 

5.15 

4.79 

4.78 

4.73 

4.72 

4.88 

Chacaltaya 

P(c)* 

13.07 

12.85 

12.75 

12.64 

12,54 

12.52 

Huancayo 

PCc)"” 

13.44 

13.24 

13.14 

13,04 

12.95 

12.92 

*No  penumbral  structure  In  the  vertical  direction  for  Chacaltaya  or  Huancayo, 
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(A)  Does  the  general  trend  of  the  vertical  cutoff  rigidity  (l.e.  Increase, 
decrease  or  no  change)  remain  constant  with  both  field  representations  for 
1980  (l.e,  using  the  original  time  derivatives  and  the  adjusted  time 
derlvat Ives) ? 

(B)  Do  significant  differences  occur  In  the  vertical  cutoff  rigidity 
values  for  the  two  field  representations?  , 


3,  RESULTS 

Table  1 presents  vertical  cutoff  rigidities  for  selected  Latin  American 
locations  where  the  vertical  cutoff  rigidity  has  been,  clecreas Ing  with  time. 

All  values  were  calculated  using  the  IGRF  Epoch  1965  field  coefficients  with 
its  associated  time  derivatives  with  the  exception  of  the  last  column,  labeled 
1980A,  where  the  adjusted  time  derivatives  were  utilized.  Table  2 presents 


TABLE  2 

VERTICAL  CUTOFF  RIGIDITIES 
FOR  SELECTED  LOCATIONS 


Cutoff  Rigidity  (GV) 


Location/ 

Epoch 

Epoch 

Epoch  Epoch 

Parameter 

1965 

1975 

1980 

1980A 

20N,  315E 

P(m) 

11.32 

11.59 

11.72 

11.72 

P(c) 

10.88 

11.45 

11.63 

11.63 

P(s) 

9.78 

10.35 

10.65 

10.64 

15N,  300E 

P(tn) 

11.08 

11.29 

11.40 

11.41 

P(c) 

10.15 

10.71 

10.97 

10.99 

P(s) 

9.10 

9.63 

9.86 

9.88 

25S,  135E 

P(m) 

8.17 

8.02 

7.91 

7.92 

P(c) 

7.93 

7.71 

7.61 

7.67 

P(s) 

6.37 

6.25 

6.35 

6.32 

30S,  330E 

P(m) 

10.29 

10.41 

10.47 

10.47 

P(c) 

7.91 

8.46 

8.64 

8.63 

P(s) 

7.49 

7.83 

7.97 

8.06 

45N,  105E 

P(m) 

6.58 

6.62 

6.63 

6.59 

P(c) 

6.27 

6.29 

6.31 

6.29 

P(s) 

5.56 

5.51 

5.18 

5.33 

25S,  195E 

P(«) 

11.69 

11.65 

11.63 

11.76 

P(c) 

11.20 

11.22 

11.23 

11.18 

P(») 

10.25 

10.29 

10.31 

10.25 

vertical  cutoff  rigidities  for  six 
locations  as  follows:  two  locations 
each  where  the  vertical  cutoff  rigidity 
Increases,  decreases,  or  remains  essen- 
tially constant  with  respect  to  time. 
The  parameters  listed  In  these  tables 
are  the  main  cone  cutoff  rigidity, 

P(m) , the  effective  cutoff  rigidity, 
P(c)  , and  the  Stormer  cutoff  rigidity, 
P(s).  The  cutoff  rigidity  values  were 
taken  from  Shea  et  al.  (1976). 


4.  DISCUSSION 

An  inspection  of  the  values  given 
In  Tables  1 and  2 show  that  the  general 
trend  (increase,  decrease  or  no  change) 
Is  preserved  for  vertical  cutoff  rigid- 
ity calculations  for  1980  irrespective 
of  the  utilization  of  the  original  or 
adjusted  time  derivatives  In  the  de- 
termination of  the  geomagnetic  field 
coefficients  for  Epoch  1980.  Dif- 
ferences In  the  actual  values  for  both 
the  main  cone  and  effective  cutoff 
rigidity  are  of  the  order  of  0.05  GV. 
(He  exclude  the  Stormer  cone  cutoff 
rigidity  In  this  discussion  since  It 
Is  the  most  difficult  value  to  accu- 
rately determine  by  the  trajectory- 
tracing  process.)  The  question  then 
arises  as  to  whether  0.05  GV  is  sig- 
nificant In  the  utilization  of  these 
values  in  analyses  of  cosmic  radiation 
data.  In  discussing  this  problem  we 
must  consider  (A)  the  relative  accuracy 
of  the  values,  (B)  the  absolute 
accuracy  of  the  values,  and  (C)  the 
uncertainties  in  the  trajectory- 
tracing  technique. 
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Relative  Accuracy.  For  analyses  where  the  cosmic  radiation  Intensity 
observations  at  one  location  are  compared  with  those  at  other  locations,  the 
vertical  cutoff  rigidity  calculated  for  1980  using  either  set  of  time 
derivatives  would  be  adequate  provided  consistency  was  maintained. 

Absolute  Accuracy.  We  know  of  no  experimental  measurements  at  the  present 
time  capable  of  determining  a 0,05  GV  difference  in  vertical  cutoff  rigidity. 

In  utilizing  vertical  cutoff  rigidity  as  an  analysis  parameter,  we  must 
remember  that  this  is  the  cutoff  value  in  one  specific  direction  (l,e.  radial 
vertical)  often  accepted  as  a value  typifying  the  location.  Absolute  dif- 
ferences of  0,05  GV  would  be  a second  order  effect  when  compared  against  the 
above  simplifying  assumption.  If  the  vertical  value  is  extrapolated  (perhaps 
by  the  application  of  Storraer  theory)  to  other  azimuth  and  zenith  angles,  then 
the  absolute  differences  of  0,05  GV  are  probably  smaller  than  the  uncertainties 
introduced  by  the  extrapolation  process. 

Computational  Technique,  The  problems  in  the  computational  technique  can 
be  separated  into  two  parts  - the  increment  of  rigidity  used  in  the  deter- 
mination of  the  vertical  cutoff  rigidities,  and  the  size  of  the  step  length 
used  in  the  trajectory-tracing  process. 

Most  vertical  cutoff  rigidities  calculated  by  the  trajectory-tracing 
procedure  have  utilized  0,01  GV  rigidity  intervals,  particularly  in  the 
penumbral  regions,  (See  Shea  et  al, , 1965  and  Shea  and  Smart,  1974,  for  a 
discussion  of  the  effects  of  interval  size,)  Even  using  a standard  interval 
of  0,01  GV  can  result  in  slight  differences  in  the  calculation  of  an  allowed 
or  a forbidden  orbit  in  the  penumbral  region  when  utilizing  two  different 
computers  simply  because  of  the  manner  in  which  each  computer  operates,  such 
as  the  number  of  significant  figures  carried,  etc.  However,  even  though  minor 
differences  in  the  determination  of  allowed  and  forbidden  orbits  in  the 
penumbra  might  be  calculated,  the  gross  features  in  the  penumbra  are  preserved. 
This  same  effect  is  true  of  the  trajectories  calculated  for  this  paper. 

Figure  1 illustrates  the  penumbral  structure  for  particles  vertically  incident 
at  15°N,  300®E  for  Epochs  1965,  1975  and  the  two  sets  of  coefficients  utilized 
for  Epoch  1980,  Although  the  effective  vertical  cutoff  rigidity  increases 
with  time  at  this  location  (See  Table  2),  the  gross  characteristics  of  the 
penumbral  sttucture  are  preserved. 

The  original  trajectory-tracing  program  of  McCracken  et  al.  (1962) 
utilized  a library  of  step  sizes  in  the  numerical  integration  technique  that 
were  deliberately  made  small  in  an  attempt  to  minimize  the  error  accumulation 
in  the  Runge-Kutta  iteration  process.  In  an  effort  to  make  the  computer 
program  more  efficient  and  less  time  consuming,  the  standard  McCracken  library 
of  step  sizes  was  replaced  by  a variable  step  size  chat  was  about  1/50  of  the 
distance  a particle  with  a specific  rigidity  traveled  during  one  gyration 
(Shea  et  al, , 1976),  Although  this  necessitated  the  recalculation  of  the  step 
length  for  each  Runge-Kutta  iteration,  the  trajectory  calculations  couid 
be  pevformed  with  considerably  greater  speed  without  loss  of  appreciable 
accuracy.  Minor  differences  between  the  two  programs  that  occurred,  primarily 
in  the  pentimbral  region,' were  attributable  to  the  slightly  different  orbits 
calculated  for  two  identical  particles. 

An  example  of  the  differences  that  may  occur  in  the  utilization  of  dif- 
ferent step  sizes  in  the  trajectory-tracing  process  is  given  in  Figure  2 which 
Illustrates  the  penumbral  structure  for  particles  vertically  incident  at 
45®S,  240“e  as  calculated  using  four  different  step  sizes  in  tracing  the  ‘ 
particle  trajectories  through  the  IGRF  (Epoch  1965)  geomagnetic  field  model. 
Although  minor  differences  in  the  fine  line  structure  can  be  ascertained,  the 
gross  features  are  similar.  These  minor  differences  result  In  the  determin- 
ation of  effective  vertical  cutoff  rigidities  that  differ  by  0,06  GV  - 
approximately  the  same  differences  that  are  found  in  the  vertical  cutoff 
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flp.iire  I,  Penumbral  structure  for 
particles  vertically  Incident  at  15®N, 
tOOOE  for  Epochs  1965,  1975,  and  the 
two  sets  of  coefficients  utilized  for 
Epoch  1980.  The  column  labeled  1980A 
indicates  Epoch  1980  utilizing  the 
adjusted  time  derivatives.  White 
indicates  allowed  rigidities  and 
black  Indicates  forbidden  rigidities. 


Figure  2.  Penumbral  structtire  for  par- 
ticles vertically  Incident  at  65**S, 
240°E,  as  calculated  using  four  differ- 
ent step  sizes;  the  original  McCracken 
step  size,  the  new  variable  step  size 
(v.  step),  and  the  variable  step  size 
divided  by  2 (v.  step/2)  and  divided 
by  4 (v.  step/4).  The  main  cone  cutoff 
rigidity  is  given  by  the  effective 
cutoff  rigidity  by  P^,  and  the  Stonser 
cutoff  rigidity  by  P^.  White  Indicates 
allowod  rigidities  and  black  Indicates 
forbidden  rigidities. 


rigidity  calculations  determined  utilizing  the  two  seta  of  coefficients  for 
a I960  Epoch  of  the  geomagnetic  field. 


5.  CONCLUSIONS 

Erom  the  results  presented  In  this  paper  we  conclude  that  the  utilization 
of  the  adjusted  time  derivatives  with  the  1975  IGRF  coefficients  does  not 
result  in  a discontinuity  In  the  calculation  of  vertical  cutoff  rigidities. 
The  cutoff  rigidity  differences  we  calculate  utilizing  the  original  time 
derivatives  for  the  IGRF  (Epoch  1965)  field  coefficients  and  the  adjusted 
time  derivatives  for  the  IGRF  (Epoch  1975)  field  coefficients  are  approximate 
ly  0,05  GV  - the  same  order  of  uncertainty  that  Is  present  In  the  computation 
al  technique  Itself.  We,  therefore,  recommend  that  the  new  adjusted  time 
derivatives  with  the  1975  Epoch  of  the  geomagnetic  field  be  utilized  In 
determining  cutoff  rigidity  values  for  the  Interval  1976-1980. 
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18  supplementary  notes 


Reprint  from  15th  International  Cosmic  Ray  Conference, 
Conference  Papers.  A,  204,  1977. 
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^0  abstract  (Conttnue  on  revarae  alda  tf  neresamry  mnd  Idenl/ty  by  block  nutuber) 

Recent  geomagnetic  measurements  have  shown  that  chanv.cs  In  the  geomagnetic  flelci 
have  deviated  from  the  predictions  of  the  original  International  Geomagnetic 
Reference  Field  (IGRF)  coefficients  and  associated  Lime  derivatives.  Accordingly 
a new  set  of  time  derivatives  has  been  determined  with  the  suggestion  that 
these  be  applied  for  time  periods  after  1975.  To  ascertain  the  effects  of  the 
changes  of  geomagnetic  field  time  derivatives  on  cosmic  ray  analyses,  vertical 
cutoff  rigidities  were  calculated  for  selected  locations  using  the  IGRF  1965 
field  coefficients  and  associated  time  derivatives  aiul  the  IGRF  1975  field 
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coefficlonts  .'md  tlio  lU'W  tlmr  <lcr  lv.it  Ivor  m.nklng  p;ic-li  sot  of  <-.-i  1 ciil  at  I ons  ap- 
proprlnto  for  a 1^)80  Kpoch.  A coinp.ir  1 son  of  those  two  sots  o'"  calculations 
shows  that  the  uncertainty  in  prcullctc'd  magnetic  fields  does  not  appear  to  be  a 
serious  limiting  factor  in  tlie  use  of  calculated  vertical  cutoff  rigidities  for 
the  analyses  of  cosmic  radiation  data. 
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